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ABSTRACT 

Detailed analysis of the substructure of Lya nebulae can put important constraints on the physical 
mechanisms at work and the properties of galaxies forming within them. Using high resolution HST 
imaging of a Lya nebula at 2: w 2.656, we have taken a census of the compact galaxies in the vicinity, 
used optical/near-infrared colors to select system members, and put constraints on the morphology 
of the spatially-extended emission. The system is characterized by (a) a population of compact, low 
luminosity (~ 0.1 L*) sources — 17 primarily young, small [Re « 1 — 2 kpc), disky galaxies including 
an obscured AGN — that are all substantially offset (>20 kpc) from the line-emitting nebula; (b) the 
lack of a central galaxy at or near the peak of the Lya emission; and (c) several nearly coincident, 
spatially extended emission components — Lya, Hen, and UV continuum — that are extremely 
smooth. These morphological findings are difficult to reconcile with theoretical models that invoke 
outflows, cold flows, or resonant scattering, suggesting that while all of these physical phenomena 
may be occurring, they are not sufhcient to explain the powering and large extent of Lya nebulae. In 
addition, although the compact galaxies within the system are irrelevant as power sources, the region 
is significantly overdense relative to the field galaxy population (by at least a factor of 4). These 
observations provide the first estimate of the luminosity function of galaxies within an individual Lya 
nebula system, and suggest that large Lya nebulae may be the seeds of galaxy groups or low- mass 
clusters. 

Subject headings: galaxies: evolution — galaxies: formation — galaxies: high-redshift 



1. INTRODUCTION 

Giant (^100 kpc), radio-quiet Lya nebulae (or Lya 
"blobs") that have been discovered in the distant uni- 
verse by virtue of their extremely luminous Lya emis- 
sion (^ 10^'^erg s~^) are thought to be regions of ongoing 
massive galaxy formation. When studied in detail, these 
systems show complex morphologies, obscured active 
galactic nuclei (AGN) and / or associat ed star-forming 
galaxies (e.g.. [Francis et al.l [1996; Ivison et al.l [T998l:_ 
iSteidel et aLll2000t lPalu~ et al. r i2004l: IChapinaiT et al.l 
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2008 ;[ 



Matsuda et g1 [20M . \2mft iBasu-Zvch fc Scharj 



Dev et all 120 05: Gc ach et al.l 120071: [Smith et al 
'rescott et aLll2009l : lOuchi et al.ll2009l: lYang et a! 



20111 ) . There i s strong evidence that the largest Lya neb- 



ulae are rare (Saito et al. 2006; Yan g et al.ll2009L 120101 : 
iPrescottI [2009 : Matsuda et al. 2011.) and typically re- 
side in the most ove r dense regions of the Univers e (e.g. 



Palunas et al! [200l IMatsuda et al.' ^200 1 [20051 [20091: 
Saito et al .1 120061 iPrescott" et al JI2OO8; Ya ng et al.l 



Unlike Lya halos observed around quasars and radio 
galaxies (e.g., [M cCarthv 1993, and refer e nces therein; 
'Wei dinger etal . 2005; iMilev et al.l 120061: iBarrio et all 
120081: ISmith et al., ,2009ir the dominant power source 
responsible for these radio-quiet Lya nebulae has 
been difficult to determine. Studies have investi- 
gated whether Lya nebulae c ould be powered by galac- 
tic superwin d outflo ws fe.g.. iTaniguchi fc Shioval [20001 : 
ITaniguchi et al. 200l lMori et al.ll2004D , photoionization 
by ob s cured AGN or star for i nation ( e. g.. [Chapman et al.l 
200I IBasu-Zvch fc Scharfl [200I IGeach et al.l 120071 . 



200a IPrescott et al.l l2009h . or gravitational coolin; 
within cold filaments ( "cold fl ows"; e.g., 'Nilsson ct al 
2OOI Smith fc Jiirvis 20071: iDiikstra fc Loe b 20Q_, 
Goerdt et all l2oToriFaucher-Giguere et all l2010f ). but 
have come to a range of conclusions. It has also been 
suggested that the large extent of Lya nebulae could 
be due to resonant scattering of Lya photons from a 
central source, with recent work providing observational 
evidence for this e ffect around Lyman-break galaxies 
(jSteidcl et al]|201l[ ). Thus, despite considerable study, 
the mechanisms responsible for the copious Lya emis- 
sion and the large extent of Lya nebulae have remained 
controversial. 
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The question of the substructure of Lya nebulae has 
received much less attention, but it is a topic that can 
provide much-needed complementary clues to the origin 
of Lya nebulae as well as to what they ultimately evolve 
into. As the potential physical mechanisms responsible 
for powering the Lya emission each have morphological 
implications, studying the morphology of the spatially- 
extended emission on kiloparsec and sub-kiloparsec scales 
can provide insight into the processes at work and the 
underlying powering mechanism in Lya nebulae. In ad- 
dition, taking a complete census of the diffuse and com- 
pact sources within Lya nebulae and studying the rela- 
tive positions and properties (luminosities, colors, mor- 
phologies, sizes) of the galaxies that reside within them is 
valuable for establishing the evolutionary state of these 
dynamic systems. Doing all of this, however, requires 
determining the membership of these crowded regions 
either with spectroscopy (typically feasible only for the 
brightest knots) or deep, high resolution imaging that 
can resolve and put constraints on the spectral energy 
distribution (SED) of faint individual sources. As yet, 
very few radio-quiet Lya nebulae have been imaged with 
HST, and those that have typically lack the critical high 
resolution constraints above the Balmer break. 

In this paper we study the sub-kiloparsec structure 
of a giant Lya nebulae at z w 2.656 using high- 
resolution imaging from HST. This Lya nebula was dis- 
covered thanks to its extreme Spitzer/MlPS 24/im emis- 
sion and its extended morphology in broad-band Bw 
imaging (LABd05: .Dev et al. 2005i . hereinafter Paper I). 
Roughly 20" (~160 kpc) in size with a Lya luminosity of 
«1.7x 10'*^ erg s~^, LABdOS rivals other known Lya neb- 
ulae in energetics and complexity. The data presented in 
Paper I revealed at least three important components 
(and potential sources of ionization) in the system: (1) 
the strong 24/im source, likely dominated by an obscured 
AGN, (2) a Lyman break galaxy to the northeast of the 
nebula, and (3) a source that does not have a counterpart 
in the ground-based imaging but that was identified near 
the center of the Lya emission due to the presence of nar- 
row, spatially unresolved IIeiiA1640 and CivA1548,1550 
emission lines in the ground-based spectrum. Follow- 
up narrow-band imaging of the surrounding environment 
revealed that LABd05 resides within a very large fil- 
amentary structure at least 50 comoving Mpc in size 
(jPrescott et al.ll2008[ ). and imaging polarization obser- 
vations have demonstrated that th e Lya emission is no t 
strongly polarized (P < 9%, 3cr: iPrescott et al.ll201l[ ). 
hinting that scattering may not be significant in this 
source. 

While revealing, our previous studies were limited 
by ground-based resolution and depth, and our result- 
ing knowledge of the system was incomplete. First, at 
ground-based resolution, it was unclear whether the sys- 
tem hosted other compact galaxies and whether the Lya 
itself contained a compact central source or was clumpy 
on small scales. Second, it was unclear whether the 
24/im source and the Lyman break galaxy were impor- 
tant power sources for the Lya nebula. The geometry 
of the system, with both the 24/im source and the Ly- 
man break galaxy offset from the centroid of the Ly a by 
2'.'5 (>20 kpc in projection), argued against this possi- 
bility, and the observed SED of both sources suggested 



that, barring inhomogeneous obscuration, they were un- 
likely to power more than ~20% of the Lya emission. 
Finally, the source of the unresolved Hen and CiV emis- 
sion, which appeared to be centered within the Lya emis- 
sion, was uncertain. While Hen and Civ emission often 
indicate shock excitation or a hard ionization source, no 
central galaxy (that could be driving shock-heating via a 
superwind) was visible in the ground-based imaging and 
furthermore the measured line ratios were inconsistent 
with shocks. 

In the present work, we use high resolution HST/ACS 
and NICMOS imaging to take a census of the compact 
sources within the system, measure their luminosities, 
morphologies, and locations relative to the line-emitting 
gas, investigate the question of the location and mor- 
phology of the Hen-emitting region, and determine the 
morphology of the Lya nebula itself. In Section [2] we de- 
scribe our observations and reductions, and in Section [3] 
we present our results on the different components of the 
Lya nebula system. Section |4] summarizes what we have 
learned about the small-scale morphology of LABdOS and 
explores the implications of these findings for our under- 
standing of what causes the Lya nebula phenomenon. 
We conclude in Section [S] In this paper, we assume the 
standard ACDM cosmology (rJM=0.3, rJA=0.7, h^O.7); 
the angular scale at z = 2.65 6 is 7.96 k pc/" . All magni- 
tudes are in the AB system (|Okelll974[ ). 

2. OBSERVATIONS & REDUCTIONS 

Paper I examined the properties of LABdOS us- 
ing Keck/LRIS spectroscopy, optical imaging from 
the NOAO Deep Wide-Field Survey (NDWFS; 
IJannuzi fc De^ IT^m. and Spitz e r/IRAC and MIPS 
imaging (EiscnhardLefaD l2004t iHouck et al] |200S). 
The large-scale environment of LABdOS was studied 
using Subaru /SuprimeCam interm ediate-band /A44S 
Lya imaging (IPrescott et al.|[2008h . and constraints on 
the Lya polarization of the system were obtained using 
imaging polarimetry with the B ok Telescope and th e 
SPOL CCD Spectropolarimeter (|Prescott et all l2011|). 
In this paper, we add high resolution HST/ACS and 
NICMOS imaging to study the small-scale morphology 
and local environment of LABdOS. Table [T] lists the 
instruments, filters, and total exposure times for the 
HST imaging. Figure [T] shows selected postage stamps 
from the multi-wavelength dataset used in this work. 

2.1. HST ACS Data 

We obtained HST Advanced Camera for Surveys 
(ACS) imaging of LABdOS on UT 2006 January 13, 
14, and 24 using the F606W^ (Veoe) filter and two 2% 
ramp filters, FRA62N (centered on Lya at zRi2.7) and 
FR601N (centered on HeiiA1640 at zw2.7)E!l Basic im- 
age calibrations (overscan, bias, and dark subtraction, 
flat-fielding) were provided by the standard HST ACS 
pipeline with On-The-Fly-Reprocessing (OTFR) and the 
task calacs. We removed a residual offset in the bias 
level of the individual amplifiers on each of the ACS de- 
tectors (roughly a 2% effect relative to the background) 
by estimating the sky background in each amplifier sep- 
arately using a sigma-clipped mean and subtracting it 

" HST Cycle 14; GO 10591 
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from the calibrated, flat-fielded individual exposures (the 
"_FLT" files). Using MuhiDrizzle's default settings and 
no additional sky subtraction, we performed the distor- 
tion correction, cosmic-ray rejection, and image combi- 
nation, yielding a final scale of 0.05"/pixel and a field- 
of-view of 207"x205". The point-spread-function (PSF) 
size for the ACS imaging is FWHM=0.07", as measured 
using the TinyTim PSF emulator The 5a point source 
limiting magnitudes for the ACS imaging are 28.3, 25.7, 
and 25.7 mag (0.4" diameter aperture) for the F606W, 
FR462N, and F601N fihers, respectively. 

The narrow-band Lya and Hen imaging contain both 
line and continuum emission. To generate HST line- 
only Lya and Hen images, we scaled the Veoe image 
by factors of 42.9 and 22.9, respectively (estimated em- 
pirically based on sources common to both images) and 
subtracted the result from the original narrow-band im- 
ages. The resulting HST line-only Lya and Hen images 
are discussed in Sections [3.2.1l and l3.2.2l To quantify the 
degree to which the Hen and Civ emission lines in turn 
contaminate the Veoe-band, we measured the total Vgoe 
flux within the same spectroscopic aperture (4.5" x 1.5") 
used in Paper I to measure the Hen and Civ fluxes, af- 
ter convolving the Vgoe image to match the ground-based 
seeing (FWHM=1"). Comparing the total Vgoe flux in 
the aperture (1.15 ± 0.01 x 10^^^ erg s~^ cm~^) to the 
reported Hen and Civ line fluxes (4.07±0.04x 10"^'' and 
4.17 ±0.04 X 10~^^ erg s~^ cm~^, respectively; Paper I), 
we find that the Heii and CiV emission lines together 
contribute a small fraction of the Vgoe-band flux within 
this aperture (< 7%), and therefore we do not apply a 
correction to the Vgoe image. 

2.2. HST NICMOS Data 

Using the NICMOS NIC2 camera on HST, we obtained 
high- resolution imaging of the source in the FllOW^ 
( Jiio) and i^ieOW^ (^leo) fllters — filters which at zw2.7 
bracket the Balmer/4000A break. The observations were 
taken during UT 2006 March 25 and 31, using a NIC- 
SPIRAL-DITH spiral dither pattern (3 point pattern 
with 0.6375" point spacing) H The da ta were reduced 
primarily using NICRED ( Magee et al.ll2007[) . Once the 
data were calibrated and corrected for electronic ghosts, 
pedestal, cosmic ray persistence, and count-rate non- 
linearity, we made the final imag e mosaics using Mul- 
tiDrizzle ([Jedrzeiewski et 311120051) and a set of custom 
bad pixel masks. The final images were supersampled to 
match the ACS pixel scale (0.05"/pixel) and have a field- 
of-view of 20"x20". The PSF sizes are FWHM=0.09" 
and 0.13", as measured using the TinyTim PSF emula- 
tor, and the 5cr point source limiting magnitudes are 27.2 
and 27.1 mag (0.4" diameter aperture) for the Jno and 
Hieo imaging, respectively. 

2.3. Subaru Suprime-Cam Data 

In this work, we make use of dee p Subaru Lya i n raging 
that was obtained previously bv iPrescott et al.l ()2008f ) 
using the Su baru telescope and th e SuprimeCam wide- 
field imager (|Miyazaki et al.l [20021 ). These observations 
used an intermediate-band filter lAMb (Ac «4458A, 

^'^ TinyTim: http:/ /www.stsci. edu/hst/observatory/focus/TinyTim"| 
15 HST Cycle 14; GO 10591 



AXpwHM ~20lA), centered on the Lya line at the red- 
shift of the nebula. The limiting magnitude of the Lya 
image is 26.6 mag (5cr in a 2" diameter aperture). Ad- 
ditional details on the observa t ions a nd data reduction 
can be found in'P rescott et al.l (j2008[) . 

We generated a Subaru line-only Lya image by sub- 
tracting off a smoothed version of the ACS Vgoe-band 
image. From the ground-based spectroscopy, we know 
that the source labeled "A" in the NDWFS Bw image 
(Figure [T]) is a Lyman Break galaxy (LEG) at the red- 
shift of the system and shows little if any Lya emission 
or absorption (Paper I) . A correct continuum subtraction 
should therefore leave the LEG with zero flux in the Sub- 
aru line-only Lya image. We smoothed the Vgoe image to 
match the PSF of the Lya image (FWHM=0.7"), resam- 
pled to the same pixel scale as the Lya image, measured 
the flux of the LEG in both the Lya and Veoe-band im- 
age (using a 1.0" diameter aperture), scaled the Veoe to 
match the flux of the LEG in the Lya image, and sub- 
tracted the two to create a continuum-subtracted Lya 
image. Since the Veoe is a rather crude approximation 
to the continuum in the Lya image, the accuracy of this 
continuum subtraction will vary with source color. We 
note that for all the galaxies at the redshift of the sys- 
tem, the subtraction should be relatively accurate; those 
that show residual emission in the Lya image are likely 
Lya-emitters themselves. On the other hand, in the case 
of a known interloper galaxy at z « 3.2 (labeled 'E' in 
Figure [TJ Paper I), we expect our continuum subtrac- 
tion procedure to overestimate the continuum (since at 
this redshift the /A445-band is sampling the continuum 
shortward of Lya which is typically depressed due to the 
Lya forest). This is consistent with the slight evidence 
of oversubtraction that we see at the position of the in- 
terloper galaxy (Figure [2|) . 

2.4. Image Registration 

To ensure accurate image registration, we generated 
catalogs of source positions in each image using SExtrac- 
tor (Bcrtin Ik Arnouts 1996) and registered all images to 
the NDWFS Bw frame using the IRAF tasks ccmap and 
ccsetwcs. We carried out the registration in three steps. 
First, the ACS Vqoq was convolved with a Gaussian kernel 
to match the PSF of the Bw image (FWHM=0.98"). We 
registered the smoothed ACS Veoe image to the NDWFS 
Bw image, applied the solution to the unsmoothed ACS 
Veoe image, and then registered the ACS Hen, NICMOS 
Jiio, and NICMOS Hiqq images to the unsmoothed ACS 
Veoe image. Finally, we registered the ACS Lya image to 
the ACS Hen image and the Subaru 7^445 image to the 
NDWFS Bw image. This sequential procedure was used 
to maximize the number of common sources available to 
compute the astrometric solution for each image pair and 
to avoid compounding registration errors. Table [2] details 
the number of sources used in the registration and the 
final estimated astrometric uncertainty relative to the 
NDWFS Bw image. The NDWFS astrometry is tied to 
a frame defined by stars in the USNO-A2.0 catalog. 

3. THE COMPONENTS OF THE NEBULA 

Our multi- wavelength observations show that LAEd05 
contains a number of compact galaxies, diffuse rest- 
frame UV continuum emission, smooth Lya emission, 
spatially-extended Hen emission, and an obscured AGN 
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- all within a «10" region. In this section, we explore 
each component of the Lya nebula system in detail. In 
Section |4l we summarize the key morphological charac- 
teristics of this system and discuss the implications of 
these results for our understanding of the physical mech- 
anisms at work in Lya nebulae, their evolutionary state, 
and the properties of the galaxies forming within them. 

3.1. Compact Sources 

To determine the relative positions of all the sources 
in the vicinity of the Lya nebula we created a composite 
stack of the Veoe, JiiO) and iJieo images, after convolv- 
ing the Vgoe and Jno-band images to the same PSF as 
the HiQo image and dividing each image by the variance 
of the sky. We generated an initial list of source posi- 
tions in the image stack using SExtra ctor (3a threshold, 
minim um contiguous area of 4 pixels; IBertin fc Arnoiitsl 
[1991 . AU sources within 7" that are detected above 
the 5(7 limiting magnitude in the Vgoe band are labeled 
with ID numbers in Figure [5] We then generated an 
additional catalog of sources positions using the uncon- 
volved Vgoe-band image and the same parameters. This 
"Vgoe-only" catalog is necessary for our analysis of the 
number counts in the vicinity of LABd05 and the associ- 
ated completeness corrections (Sections 13. 1.11 and 14. 2. ip . 
In four cases (#1, 21, 26, and 46) where SExtractor did 
not deblend an apparent close object pair, we manually 
added a second source position to the catalog (#57, 58, 
59, and 60, respectively), as described in Section [3.1.41 
The current data are not sufficient to distinguish whether 
these objects are true companions or just morphological 
peculiarities (i.e., tidal features, dust lanes, etc.) as- 
sociated with the primary object. We choose to treat 
these pairs as separate objects; however, combining them 
does not significantly change our conclusions. Individual 
postage stamps extracted from the image stack as well 
as the Veoe, >^iiOj and ffieo-band imaging are shown for 
all sources in Appendix |^ 

3.1.1. Optical and Near-Infrared Photometry 

We measured aperture photometry (0.4" diameter 
apertures) in all three bands using the original uncon- 
volved images and the positions derived from the image 
stack. The aperture size was chosen in order to contain 
as much flux as possible while minimizing contamina- 
tion from neighboring sources. Aperture corrections of 
[1.15, 1.52, 1.76] were computed using the TinyTim PSF 
emulator and applied to the Veoe, Jiio, and iJieo pho- 
tometry, respectively. The resulting photometry is given 
in Table 13 

Without knowing the intrinsic colors of sources as a 
function of magnitude, it is difficult to estimate the com- 
pleteness of the stacked catalog. Instead, we measured 
the completeness in the Vgoe-band alone using the stan- 
dard approach. First, we generated simulated galaxies, 
modeled as Gaussian profiles with Rg = 0.12" (1 kpc, 
the typical size of non-member sources in the field; Sec- 
tion 13.1.41) . We then inserted the simulated galaxies 
into the Vgoe-band imaging. We generated a "Veoe-only" 
source catalog in the same manner as above and com- 
puted the completeness of our approach as a function of 
input Vqo6 magnitude. The 80% and 50% completeness 
limits are 27.6 and 28.0 mag in the Vgoe band, respec- 
tively. 



3.1.2. The Obscured AGN 

Paper I postulated the existence of an obscured AGN 
at the position of the MIPS source based on the strong 
24/im emission and the shape of the full SED. The 
corresponding detections in the IRAC bands agreed 
with the MIPS source position to within the astro- 
metric uncertainty (w0'.'5) and showed that the source 
had a power-law SED in the mid-infrared, character- 
istic of an obscured AGN (e.g., Alonso- Herrero et al.l 
I2006D . Later IRS spectroscopy demonstrated that the 
infrared source is a t the redshift of the Lya nebula 
(jColbert et al.|[20TlT ). and millimeter and submillimeter 
observations have confirmed that the full SED is best 
appro ximated by a Mrk 231 (i.e., AGN-domina ted) tem- 
plate (|Bussmann et al.l 120091: lYang et al.ll2012D . In Pa- 
per I, the physical location of the AGN within the system 
was uncertain due to the lower resolution of the IRAC 
and MIPS imaging, but the centroid of the mid-infrared 
emission appeared to be offset to the north of the bright- 
est Lya emission. The addition of the HST/NICMOS 
imaging revealed an extremely red source (#36) located 
at the centroid of the MIPS 24/im emission (Figure [Ij 
that is very centrally concentrated (Section 13. 1.4p . This 
source shows a strong Balmer/4000A break — it is barely 
detected in Vgoe and Jno but very bright in -ffieo — and 
is one of the 17 sources fiagged as members of the system 
(Section [3.1.3p . It is located in a crowded region, with 5 
close neighbors within sal. 5" («12 projected kpc), and 
diffuse emission visible in the NICMOS i?i6o band, sug- 
gestive of an ongoing merger. Since it is plausible to 
assume that the AGN lies near the deepest part of the 
gravitational potential well of this system, we will take 
the position of source #36 as the center of the system for 
our subsequent analysis. The measured projected offset 
between the AGN and the centroid of the Lya emission 
(Section [SXI]) is «1.9" («15 kpc). 

3.1.3. Assessing System Membership 

Only two compact sources were previously identified 
from ground-based imaging (Paper I; see NDWFS Bw 
image in Figure[T]). At the Northeast corner of the system 
is a compact source — labeled "Galaxy A" in Paper I — 
that ground-based spectroscopic follow-up showed to be 
an LBG at the redshift of the system. The source at the 
Southwest corner — named "Galaxy B" in Paper I — was 
argued to be an interloping system based on the identifi- 
cation of Lya at z w 3.2 in the ground-based spectrum. 
The high resolution HST imaging resolves both of these 
objects into two components: Galaxy A is associated 
with objects #26 and 59 and Galaxy B contains knots 
#46 and 60 (see Section 13.1.41 and Appendix . Since 
the spectroscopic identification was done using ground- 
based spectroscopy that was unable to resolve the two 
components in each case, it is possible that these pairs 
are in fact due to chance coincidence. However, with sep- 
arations of only ~ 0.2" (i.e., 1.6 kpc), chance projection 
is extremely unlikely, and the photometry shows that the 
colors of both components in each pair are similar. Given 
the ground-based spectroscopic redshift, the very small 
likelihood of a chance coincidence, and the similar colors, 
we will assume from here on that #26 and 59 are both 
associated with the Lya nebula system at z « 2.7 and 
that #46 and 60 are both interlopers at z ~ 3.2. 



Resolving the Galaxies within a Giant Lya Nebula 



5 



To determine the membership of the remaining 
sources, we make use of the measured optical/NIR col- 
ors. Since the NICMOS Jno and -ffieo bands straddle 
the Balmer/4000A break at the redshift of the nebula, 
sources within the system should show red Jhq ~ iJigo 
colors if they have evolved enough with time. While a full 
SED-fitting approach would be poorly constrained with 
only three bands, we can use this fact to identify other 
sources that are likely associated with the system. We 
start by selecting the sample of sources that are within a 
radius of 7" from the AGN (#36) and that are brighter 
than the 5a limiting magnitude in all three bands. In 
Figure 12] we plot the Veoe — "^iio versus Jno — -ffieo col- 
ors of the resulting sample along with greyscale contours 
representing the expected color distribution of galaxies 
drawn from the field (taken from the Hubble Ultra Deep 
Field, HUDF, w hich used the sa me instrument and fil- 
ters as this work: lCoe et al.ll2006l ). The subset of HUDF 
galaxies with photometric redshifts consistent with the 
systemic redshift to within typical photometric redshift 
errors {Zphoto = 2.656 ± 0.15) is shown with line con- 
tours. While a handful of sources have colors entirely 
consistent with being drawn from the field (i.e., they are 
near the peak of the HUDF greyscale), there is also a 
locus of objects that is broadly consistent with being at 
the systemic redshift (i.e., they are within the line con- 
tours) but that extends along a line roughly parallel to 
the reddening vector towards much redder colors. These 
red colors are quite unusual for typical field galaxies, a 
hint that this locus may be composed primarily of system 
members with varying amounts of dust. 

In Figure m we plot the same sample alongside a series 
of age tracks for simple stellar population models (single 
unreddened bursts, solar metallicity) at z = 1.5 — 4.0 as 
well as a constant star-forming model (solar metallicity) 
and a low metalli city model (Z = 0.0001), both at the 
systemic redshift ()Bruzual fc Chariot! I2003D . We again 
see that the locus discussed above consists of sources 
that are consistent with the systemic redshift if we al- 
low for a low to moderate amount of dust extinction 
{E{B — F) « 0.0 — 0.4 mag). We draw a dividing hue in 
color-color space with a slope parallel to the reddening 
vector in order to select sources that are consistent with 
the systemic redshift, given typical photometric errors. 
Those to the upper left of the line we consider "members" 
and those to the lower right, "non-members" (designated 
"NM" ) . As a check on the effectiveness of our approach, 
we apply the same color cut to the HUDF galaxy cata- 
log and plot a histogram of the photometric redshifts for 
this sub-sample (Figure [S]). The color cut is effective at 
selecting high redshift galaxies (90% are at Zphoto ^ 1), 
and in contrast to the redshift distribution for the full 
HUDF sample, the photometric redshifts of HUDF galax- 
ies selected using this simple color cut are peaked at the 
systemic redshift of LABd05 {zphoto ~ 2.66). On the 
other hand, our color cut is clearly approximate. The 
comparison sample selected from the HUDF contains a 
small subset (10%) of galaxies with Zphoto ^ 1, which is 
consistent with the fact that for young (w 25 Myr) sin- 
gle burst models the predicted colors at z < 1 overlap 
those for higher redshifts. The color cut is also not able 
to reject the known interloper system a.t z 3.2 (#46 
and #60, or Galaxy B). From Figure [Sj this is not at all 
surprising, as the peak of the Zphoto distribution is broad, 



spanning 2.1 < Zphoto ^ 3.2. Furthermore, the fact that 
the interloper system is located at the young end of the 
age tracks in color-color space (Figure |4|) is consistent 
with the detection of Lya emission (Paper I). 

We assess the robustness of our membership assign- 
ment further in Section r4.2.2l At this point, we subdivide 
the member sample based on the projected distance from 
the AGN (#36). The 9 sources within the radius of 2.1" 
(chosen to include all spectroscopically confirmed mem- 
bers, #26, 59, 36) are considered to be members with 
high confidence and designated "Ml". The 8 members 
outside this radius are designated "M2" . In Figure El we 
see that all the member sources (circled) are offset by 
1.9 — 6.7" (15 — 53 projected kpc) from the peak of the 
Lya emission (Section I3.2.ip . i.e., they effectively reside 
at the outskirts of the Lya nebula. 

3.1.4. Sizes and Morphologies 

Our next step was to use GALFIT (jPeng et al.ll2002[ ) to 
derive sizes and morphologies for the sources within the 
vicinity of LABd05. Visual inspection of the Vqoq image 
revealed that in addition to the many compact galaxies in 
the region, there is diffuse, spatially extended emission. 
In order to avoid biasing the fits for the compact sources, 
we began by subtracting off an approximate fit to this 
diffuse emission using GALFIT. We then fit all remain- 
ing compact sources within 7" of the AGN that have a 
peak surface brightness brighter than 25.3 mag arcsec"^ 
and an isophotal magnitude brighter than 29.0 AB mag 
in the T^goe band. To speed up the process of fitting so 
many sources simultaneously, we performed the fits in 
small batches initially, i.e., fitting 1 — 7 nearby sources 
at a time while masking the remaining sources. After 
experimentation showed that sources #1, 21, 26, and 46 
were not being fit well using single components, sources 
#57, 58, 59, and 60, respectively, were manually added 
to the catalog. During the initial batch-fitting, the posi- 
tion of sources #57, 58, and 59 were then masked, while 
#60 was fit normally. Once appropriate fitting parame- 
ters had been determined for all sources, we performed 
a final run of GALFIT on the original image fixing all 
parameters for the compact sources to the pre-derived 
values. We allowed the few remaining special cases (#57, 
58, and 59) as well as the diffuse continuum component 
to be fit freely by GALFIT. The entire fitting process 
was repeated using the fitting parameters derived for the 
diffuse continuum component in this final step as the ini- 
tial guess for the first step. The resultant composite fit 
and corresponding residual images are shown in Figure [71 
The centroids and morphological results for the compact 
sources derived from the Vgoe band are given in Table |3l 
and for the diffuse component in Table [H We discuss 
this diffuse UV continuum component in further detail 
in Section [3.2.31 

3.2. Diffuse Components 

In addition to the population of compact sources, 
LABd05 hosts several nearly coincident diffuse emission 
components. 

3.2.1. Diffuse Lya Emission 

If the Lya nebula is powered by sources within the 
nebula itself, we would expect that to be reflected in the 
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morphology of the Lya emission. We used the high res- 
olution ACS Lya imaging to look for spatially resolved 
knots or clumps that could signal the locations of the ion- 
izing sources for the nebula. However, we find that the 
morphology of the Lya nebula is smooth, with no sig- 
nificant substructure (Figure [8|) . At the depth of these 
observations we would have detected point-source regions 
down to a Lya flux of F^ya = 2.6 x 10~^^ erg cm~^ 
(5cr in 0.4" diameter apertures) or Lya luminosity of 
LLya = 1-5 X 10^^ erg s~^. This corresponds to a 
star-formation rate of 1.3 Mp^ yr~^ (assuming Case 
B, LLya = 8.7 X LHa', lKennicut"g I1998D . We can also 
rule out high surface brightness clumps within the cloud 
down to a peak surface brightness limit of «4.0xl0~^^ 
erg s~^ cm~^ arcsec~^(3CT). 

Deeper Lya imaging at ground-based resolution ob- 
tained with the Subaru Telescope showed that the Lya 
emission on larger scales has a smooth elliptical morphol- 
ogy that is well fit by an exponential disk profile. Using 
GALFIT, we derived a centroid, effective radius, Sersic 
index, position angle, and axis ratio for the Lya emis- 
sion (Table SI . The resulting surface brightness profile of 
the Lya emission is shown in Figures l9l and [TOl Surface 
brightness profiles were also computed for both the Sub- 
aru Lya imaging and the ACS Lya imaging using the 
IRAF task ellipse, where the ellipse parameters were con- 
strained to the centroid, single position angle, and axis 
ratio of the GALFIT parametric fit to the Subaru Lya 
imaging. In addition, we measured aperture photome- 
try at the center of the nebula using a small aperture 
(1.6" diameter) chosen to minimize contamination from 
neighboring compact sources (Table |4]). 

3.2.2. Dijfuse He II Emission 

From ground-based spectroscopy, LABdOS is known to 
have strong Hen and CiV emission (Paper I) located near 
the center of the Lya emission ([-0.9", -2. 5"]; Figure [TT]). 
The goal of our Hen imaging was to better localize the 
Hen-emitting region within the system, and given the 
depth of our observations, we would have detected the 
Hen emission at a SNR w 10 if it were emitted as a point 
source. However, the ACS Hen imaging showed no detec- 
tion down to a 5a point source limiting flux of 1.9 x 10~^^ 
erg s^^ cm^^ (0.4" diameter aperture) and a Scr surface 
brightness limit of 3.3 x 10~^^ erg s~^ cm~^ arcsec"^ 
(Figure [TT|). While we cannot, therefore, pinpoint the lo- 
cation of the Hen emission beyond what was determined 
in Paper I, this non-detection does put a constraint on 
the size of the Hen-emitting region. To quantify this, we 
inserted a series of simulated Hen sources, modeled as 
Gaussian profiles scaled to match the measured Hen flux 
and a range of FWHM sizes of 0.1 — 1.0", into the Hen 
image (Figure [TT]) . For each of 100 Monte Carlo trials, 
we measured the observed flux of the simulated source 
) as a function of FWHM and determined the 
FWHM for which the simulated source is detected at 
F simulated I sky = 3, whcrc asky IS thc Icr limiting flux 
of the Hen image (0.4" diameter aperture). We con- 
clude that the source of Hen must be extended with a 
FWHM > 0.58" in order to be undetected in our Hen 
imaging. At the same time, the fact that the previous 
long-slit spectroscopy did not resolve the Hen line puts 
an upper limit of ~1" on its true spatial extent (Pa- 



per I). This size range corresponds to a Hen-emitting 
region that spans ~ 4.6 — 8.0 kpc at the redshift of the 
nebula. These constraints are summarized in Tabled 

3.2.3. Diffuse Continuum Emission 

The broad-band ACS Vgoe data revealed diffuse emis- 
sion located near the center of the Lya nebula (Figure [71 
Section l5.1.4[) . As Hen and Civ emission at the systemic 
redshift are located within the Vgoe bandpass, and as the 
band also straddles Lya at the redshift of Galaxy B (the 
interloper system at z « 3.2), we must flrst address the 
question of whether this diffuse Veoe-band emission could 
be dominated by line emission. However, as we showed 
in Section I^TTl the Hen and Civ emission lines can con- 
tribute at most 7% of the Vgoe-band flux. The Lya emis- 
sion from Galaxy B {FLya — 5.15 x 10~^^ erg s~^ cm~^; 
Paper I) can contribute at most another 4% of the Veoe- 
band flux, but even this negligible fraction is likely a 
gross overestimate. The spectroscopic data shows that 
the Lya emission from Galaxy B is compact, centered 
on the interloper system, rather than spatially extended 
(Paper I). We conclude therefore that the diffuse Vgoe- 
band emission seen in LABd05 is indeed rest-frame UV 
continuum emission. 

As described in Section 13.1.41 we measured the cen- 
troid, size, and luminosity of this component using 
GALFIT, avoiding contamination from nearby compact 
sources during the fitting process by constraining all 
the fitting parameters for the compact sources to the 
previously-derived values. Table [4| gives the position, 
magnitude, effective radius, Sersic index, position an- 
gle, and axis ratio of the diffuse component. The dif- 
fuse continuum component is well-fit by an exponential 
disk profile and the centroid is nearly coincident with 
that of the Lya emission (offset by « 0.8"). We mea- 
sured aperture photometry at the center of the nebula 
in all three bands using a small aperture (1.6" diame- 
ter) chosen to minimize contamination from neighboring 
compact sources (Table[4]). The surface brightness profile 
of the diffuse continuum emission was also computed us- 
ing the IRAF task ellipse, where the ellipse parameters 
were constrained to the centroid, single position angle, 
and axis ratio of the GALFIT parametric fit (Figures [Hi 
and[ini). 

At the depth of these observations we would have 
detected point-source regions with fiuxes of 3.5xl0~^^ 
erg s~^ cm~^ (5(7, 0.4" diameter aperture), correspond- 
ing to a stellar mass limit of w 3.1 x 10^ Mp^, assuming 
a 25 Myr simple stellar population (iBruzual fc CharlotI 
[200l . We can also rule out high surface brightness 
clumps down to a peak surface brightness limit of 
«9.5xl0~^^ erg s~^ cm~^ arcsec~^(3tT). However, it is 
worth remembering that if there are many lower mass 
sources distributed over this area, they would appear as 
unresolved diffuse continuum emission. 

4. DISCUSSION 

4.1. The Substructure of LABdOS 

The high resolution imaging of HST provides a pre- 
cise look at the sub-kiloparsec morphology of LABdOS, 
with important implications for our understanding of the 
substructure of Lya nebulae and the properties of the 
galaxies forming within them, as well as for determining 



Resolving the Galaxies within a Giant Lya Nebula 



7 



what power sources are ultimately responsible for the 
Lya emission. The key morphological characteristics of 
this Lya nebula system are the following: 

• Many Compact, Low Luminosity Galaxies: The 
system hosts 17 primarily small, disky, low lumi- 
nosity galaxies including an obscured AGN. 

• Offset Morphology: All the compact sources within 
the system are located >20 kpc away from the peak 
and centroid of the Lya nebula. 

• No Central Galaxy: LABdOS has no central galaxy 
or compact source brighter than w0.03 L* visible 
within the highest surface brightness region of the 
Lya nebula. 

• Diffuse Line and Continuum Emission: While the 
compact sources appear to avoid the region of the 
nebula, there are three nearly spatially coincident 
extended emission components located near the 
center of the nebula: Lya, Hen, and rest-frame 
UV continuum emission. 

• Smooth, Non-filamentary Morphology: Both the 
Lya emission and the diffuse UV continuum emis- 
sion are smooth and show surface brightness pro- 
files that are consistent with exponential disks 
(n ~ 0.7 — 0.8) with moderate axis ratios (6/a ^ 
0.7 — 0.8), i.e., they are not particularly clumpy or 
filamentary. 

• Similar UV and Lya Surface Brightness Profiles: 
The Lya and UV continuum surface brightness 
profiles, while not identical, are comparable in 
shape and extent. 

In what follows, we discuss each of the key morpho- 
logical findings in detail. Focusing first on the compact 
sources, we show that the LABd05 system is overdense 
relative to the field, hosting a population of small, disky, 
low luminosity galaxies that, while numerous, are ir- 
relevant to the ionization of the nebula (Section 14. 2|) . 
We then compare the remaining observed morpholog- 
ical characteristics to the expectations for Lya nebu- 
lae driven by superwind outfiow, cold flow, and reso- 
nant scattering scenarios and find significant discrepan- 
cies (Sections 14. 31 - 14. 6p . We end with a discussion of the 
possibility that the LABdOS system is a forming galaxy 
group (Section l4Jl) . 

4.2. Many Compact, Low Luminosity Galaxies 

4.2.1. An Overdense Regton 

Even without determining the membership of individ- 
ual nearby galaxies, it is clear from the high resolution 
ACS and NICMOS imaging that there are a large num- 
ber of compact sources in the vicinity of LABdOS, more 
than would be expected for a region of this size in the 
field. Figure [T^ shows the Vgoe-band number counts for 
sources detected above the 5a limiting magnitude and 
located within a 7" radius of the AGN in LABdOS. For 
comparison we show number count measurements drawn 
from our entire 207" x 20S" ACS pointing as well as 
field measurements taken from the Hub ble Deep Field- 
North and -South (HDF-N, HDF-S: Wilhams et al.lll996l: 



Casertano et~al| [2000l). the HUDF (jBeckwith et al.|[200l 
ML 



Coe et al.M2006D. and the GOODS-North and -South 
fields (|Giavalisco et all I2004D . Dividing the LABdOS 



number counts by those from the HUDF shows that the 
Ly a nebula system is overdense by at least factor of 4 ± 1 
relative to the field at magnitudes of Vgoe < 27. S (Fig- 
ure ll3p . with the uncertainty computed assuming Poisson 
statistics. If we further restrict the region of interest to 
a 2.S" radius region encompassing just the group of ob- 
jects lying immediately to the north of the nebula (i.e., 
centered at [0",-t-l.S"] in Figure [5]), the overdensity is a 
factor of 1S±S above the field. We can derive a better es- 
timate of the true overdensity factor relative to the field 
at this redshift by applying the membership color cut dis- 
cussed in Section [3X3l to both the LABdOS and HUDF 
galaxy catalogs. In this case, the overdensity factor of 
the 7" radius region approaches 12 ± 3 at magnitudes 
of Veoe < 27.5. These high overdensity factors are not 
surprising if LABdOS is indeed a region of active galaxy 
formation destined to become a galaxy group or cluster. 

4.2.2. Properties of the Member Galaxies 

The previous section looked at the question of whether 
this region is overdense in a statistical sense relative to 
the field, i.e., without drawing on any knowledge of sys- 
tem membership. In Figure [TH we again show the Veoe 
number counts (top panel), but in addition, we make 
use of our membership assignment to plot the number 
counts for sources that are likely within the system (bot- 
tom panel), i.e., the observed luminosity function. We 
find that the observed luminosity function of the M1-I-M2 
subset is largely consistent with what we would predict 
statistically by taking the observed Vgoe number counts 
for the LABdOS region and simply subtracting off the ex- 
pected Vgoe number counts for the field, as derived from 
the HUDF (Figure [TH dashed line). This good agree- 
ment between the statistical treatment of the previous 
section and the individual object membership treatment 
gives us confidence both that our membership assign- 
ment based on optical/NIR colors is reasonable and that 
our estimate of the observed luminosity function of the 
system is relatively robust even if our membership as- 
signment is not perfect on an object-by-object basis. To 
our knowledge, this is the first time a luminosity function 
has been estimated for an individual Lya nebula system, 
and it implies that LABdOS is dominated by low lumi- 
nosity galaxies. Assuming that M* at z w 2.7 is « —20.8 
(AB mag at rest-frame 1700A: lReddv et al.|[2008l ). which 
corresponds to m* « 24.49 AB in apparent Vgoe magni- 
tudesF^ we find that the LABdOS system is dominated 
by ~ Q.IL* galaxies. Only one galaxy, the LBG (#26), 
is more luminous than « 0.4L*. 

These data also reveal clues about the ages and mor- 
phologies of the galaxies within LABdOS. The colors of 
sources within the system appear to be consistent with 
young ages (w 2S — 100 Myr) and a range of dust ex- 
tinctions {E{B - V) ^ 0.0 - 0.4). The distribution 
of sizes and Sersic indices of the galaxies in the vicin- 
ity of LABdOS are shown in Figure [TSl with those that 
are likely associated with LABdOS shown as filled (sub- 
set Ml) and hatched (subset M2) histograms. In gen- 

The central wavelength of Veoe is Ac ~ 5907A observed Ri 
1600 A in the rest frame at z fa 2.7 
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eral, the galaxies in the vicinity of LABdOS are small 
{Re = 0.5 — 3 kpc) with Sersic indices clustering around 
n ~ 1, indicative of exponential disk morphologies. 

In summary, our analysis suggests that the LABd05 
system hosts a large population of young, small, disky, 
low luminosity galaxies. 

4.2.3. Energy Budget of the Member Galaxies 

The presence of a large population of compact galax- 
ies within the LABdOS system raises an obvious question: 
are they responsible for powering the Lya emission via 
photoionization? We find that the total ionizing lumi- 
nosity (L ■> ) of all 16 of the well-detected nearby 

galaxies at the system redshift (excluding the obscured 
AGN, #36) is 8.5 x 10^3-6.6 x 10^^ photons s"^ (unred - 
dened, 5-25 Myr single bursts iBruzual &: Chariotll2003[ ). 
ignoring the effects of distance, geometric corrections, 
and non-unity escape fractions. This is much less than is 
required to power the observed Lya (1.7 x 10^^ photons 
s~^; Paper I). The fact that all of these galaxies lie tens 
of kiloparsecs away from the peak of the Lya emission 
makes their potential contribution even smaller. 

The only single compact source that could potentially 
contribute substantially to powering the nebula is the 
AGN (#36), but simple energetic arguments imply that 
even it may not be powerful enough to explain all of 
the Lya emission. In Paper I, we found that the total 
ionizing luminosity of the AGN, estimated from an ex- 
trapolation of the mid- infrared SED, was w 1.8 x 10^^ 
photons s"^ and potentially orders of magnitude larger 
than the contribution from the other compact sources. 
There is a large associated uncertainty on this estimate 
due to the difficulty of extrapolating the ionizing lumi- 
nosity from the mid-infrared SED, but if it is accurate, 
it implies that the AGN can contribute at most 18% of 
the necessary io nizing photons to exp lain the Lya (Pa- 
per I). Recently, iColbert et al.l ()201ll) argued, based on 
measurements of PAH features in IRS spectroscopy, that 
roughly half of the bolometric luminosity of this mid- 
infrared source may arise from star formation rather than 
AGN activity. If so, our estimate of the fraction of ioniz- 
ing photons contributed by the AGN must be reduced by 
an additional factor of 2. Thus, for the AGN to be the 
dominant power source for the Lya nebula, the geome- 
try of the system must be such that the AGN is highly 
obscured to our line-of-sight but relatively unobscured 
in the direction of the gas cloud. A full analysis of the 
energetics of the system and the dominant power source 
will be addressed in an upcoming paper. 

4.3. An Offset Morphology with No Central Galaxy 

A key observational result of this work is that there 
is no galaxy or AGN at or near the center of the Lya 
nebula itself. All the compact sources visible in the HST 
imaging are offset by w20 kpc from the peak of the Lya 
emission. In fact, the member galaxies actually appear 
to encircle the nebula (e.g., #1, 18, 21, 26, 29, 45, 52, 
57), although given the small sample size, it is not clear 
whether or not this effect is significant or merely coin- 
cidental. Neither the Lya nor the diffuse UV contin- 
uum components show any significant knots of emission 
that might represent bright star-forming galaxies or lumi- 
nous stellar clusters. The lack of compact components 



brighter than 3.5x10^^^ erg s^-^ cm^^ (5a) in the dif- 
fuse Veoe-band light implies that there are no compact 
sources with luminosities brighter than ~ 2 x 10^^ erg s""'^ 
(0.03L*) embedded within the central line-emitting re- 
gion of the nebula. Similarly, the lack of any obvious 
substructure in the ACS Lya image suggests that there 
are no compact sources forming stars at a rate greater 
than ~ 1 Mq yr~^. The lack of 24/im emission spatially 
associated with the Lya nebula suggests there is also 
no hidden starburst at the center. This lack of a cen- 
tral galaxy has important implications for the potential 
power source in LABd05. 

An early model that was proposed to explain the Lya 
nebula phenomenon was shock-heating in starburst su- 
perwinds. However, from previous work we know that 
the Civ and Heii line ratios in LABd05 are not consis- 
tent with shock ionization (Paper I). More importantly, 
with no central source present to drive a wind it is ex- 
tremely difficult to describe this system using a simple 
outflow model. An outflow driven by one of the other 
sources would need to be severely asymmetric (no sec- 
ondary Lya peak is present on the opposite side) and 
extend more than ^40 kpc in one direction. 

The lack of a central galaxy is similarly problematic 
for models that attempt to explain large Lya nebulae 
such as LABd05 as evidence for cold flows. In this 
scenario, the Lya is powered primarily via gravitational 
cooling radiation, as gas falls into a gravitational 
potential well along cold filaments, heats collision- 
ally, and cools via Lya. The theoretical expectation, 
borne out more recently in numerical simulations, is 
that the extended Lya emission surrounds a grow- 
ing galaxy jH aiman et al. 2000; Fardal ct al. 200fl 
Furlanetto et al. 2005; Yang ct al. 2006; Diikstra et ahl 
2006allbl : iGocrdt ct al. 2010; Faucher-Gigucre et al.l 
2010f ). This basic prediction of the cold flow model is 
violated in the case of LABd05. In order to reconcile the 
data with this scenario, we would have to conclude that 
we are witnessing the Lya nebula at an epoch before the 
star clusters within it have dynamically relaxed to form 
a centrally concentrated galaxy. Even in this case, since 
the observed centroid of the Lya emission in this model 
should be a direct indication of the position of the center 
of mass of the system, we would need to understand how 
all of the compact sources — the luminous, obscured 
AGN and the 16 other compact galaxies — ended up 
> 20 kpc away from the center of mass. 

Finally, another hypothesis for explaining the observed 
properties of Lya nebulae is one in which Lya photons 
produced by a central source are resonantly scattered 
in the surrounding gas out to much larger radii. Re- 
cent observations have uncovered what appear to be Lya 
scattering halos around continuum-selected LBGs, with 
Lya surface brig htness profiles simil ar in shape to those 
of Lya nebulae (jSteidel et al.l 1201 1|) . and have revived 
the long-standing question of whether Lya nebulae could 
simply be scaled up versions of this phenomenon. In ad- 
dition, imaging polarimetry of one of the z = 3.1 Lya 
nebulae in the SSA22 field (LABI) res ulted in a detectio n 
of polarization of the Lya emission (jHaves et al.ll201l[ l. 
suggesting that in this case Lya photons may be scatter- 
ing from central sourc es embedded within the nebula (cf. 
iWeiimans et al.ir2010f ) rather than being produced in situ 
at large radii. While resonant scattering must contribute 
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at some level to the overall extent of Lya nebulae, it is 
not yet clear whether this process is the dominant cause 
of the large sizes observed in all (or even most) Lya 
nebulae. Indeed, this explanation is difficult to reconcile 
with our observations of LABdOS. While existing imag- 
ing p olarimetric constraints for LABd05 (Prcsco tt et alj 
1201 ll ) are not sufficient to rule out scattering entirely, 
with no central galaxy or galaxies, there is no obvious 
source of Ly a photons that could undergo resonant scat- 
tering to produce the observed Lya nebula. Instead, the 
Ly a photons could only be supplied by the AGN or com- 
pact galaxies at the outskirts. The fact that the Lya is so 
dramatically offset from these potential source(s) of Lya 
photons would require that either the gas distribution or 
the illumination is severely asymmetric. 

We note that LABdOS is not the only Lya nebulae 
known to lack a central source. The extensive multi- 
wavelength imaging of the GOODS-S field revealed a 
Lya-emitting nebula at z w 3.16 with no obvious contin- 
uum counterparts. This fact was used by the authors to 
argue the Lya emission is r nost likely powered by grav- 
itational cooling radiation (jNilsson et al.l 120061 ) . which, 
if true, would again require the system to be in a very 
early state prior to the onset of significant star formation 
at or near the center of mass of the system. However, 
the deep HST multiband imaging of the field does reveal 
a number of compact galaxies within 1-7" of the Lya 
emission peak. While none of the available photometric 
redshifts for these galaxies closely matches the redshift of 
the Lya nebula, four of the brighter cases are consistent 
with the systemic redshift, given the large redshift error 
bars, and a number of fainter galaxies without reliable 
photometric redshifts are located within 1-3" of the Lya 
emission peak. In the context of what we have learned 
about LABdOS, it is clear that powerful Lya nebulae can 
exist substantially offset from all associated continuum 
sources. We argue therefo r e tha. t it remains to be seen 
whether the iNilsson et al.l (|2006f) Ly a nebula is in fact 
alone. 

4.4. Diffuse Line and Continuum Emission 

The ACS and ground-based broad-band images show 
clear evidence for a spatially extended diffuse UV con- 
tinuum component that is co-located with the spatially 
extended Lya emission. In addition, the non-detection 
in our narrow-band FR601N imaging strongly suggests 
that the Hen emission is spatially extended as well 
(0.58" < FWHM < 1.0"). What is the origin of these 
diffuse components? Given the lack of a central galaxy, 
there are only two plausible sources: in situ spatially dis- 
tributed star formation and/or the obscured AGN that 
lies > 20 kpc away. 

One intriguing possibility is that the diffuse UV con- 
tinuum emission is due to star-formation taking place in 
very small, widely distributed, perhaps dynamically un- 
relaxed regions that are unresolved and unseen by the 
current HST imaging. Possible evidence for this type of 
extended star formation has been seen in another radio- 
quiet Lya nebula (SSA22-LAB1: IMatsuda et al.l 120071) 
as well as in the outskirts of a radio galaxy (MRC 1138- 
262, "the Spiderweb galaxv" : iHatch et all [20081). This 
scenario would explain the morphology of the observed 
diffuse continuum component and the presense of ion- 
izing radiation emerging over an extended region. The 



roughly elliptical shape of the UV continuum and Lya 
components, along with the velocity profile of the Lya 
reported by Paper I, lead us to speculate that we could be 
observing a large inclined disk exhibiting solid-body rota- 
tion. Using the central aperture measurements for both 
components (Table S]) , we derive similar star formation 
rates from the Lya and diffuse UV continuum emission 
(8.3 ± 0.1 and 5.6 ± 0.5 Mq yr~^, r espectively, within 
the central 2 arcsec^: rKennicut^ll998[) . The correspond- 
ing rest-frame equivalent width for Lya of pa 200 ± 18A, 
which is plausible for a stellar population ([Chariot fc Falll 
[1995 [Malhotra fc Rhoadl [2001 [Schaereii 120031) . How- 
ever, the total diffuse UV continuum flux implies that 
the expected contribution to the ionizing photon flux 
is 1.7 X 10^"^ — 1.3 X 10^^ photons s^^ from a young 
(5 — 25 Myr), solar metallicity, spatially distributed stel- 
lar population with a mass of 5.0 x 10^ — 4.6 x 10 ^ 
Mq (unreddened single bursts: iBruzual &: Charlotl2003() . 
This is only a fraction of the ionizing flux required to 
power the entire Lya nebula (1.6 x 10^^ photons s~^; 
Paper I), suggesting that a much lower metallicity stel- 
lar population would be required in this scenario. The 
inferred spatial extent of the He ii emission could poten- 
tially support of this picture; however, this would imply a 
stellar population of extremely young age and low metal- 
licity (< 2 Myr, Z < IO-'^Zq; Schacrer"2003). 

The alternative is that the diffuse line and continuum 
emission are the result of photoionization and scattering 
from the AGN located nearly 20 kpc away from the peak 
of the diffuse line and continuum components. Evaluat- 
ing this scenario requires knowledge of the total power of 
the AGN and the degree of obscuration in the direction of 
the nebula. The bolometric luminosity derived from the 
infrared SEP ind icates a powerful AGN ( 8.6 x lO^^LfT,; 
lYang et all 120121 ). On the other hand, IColbert et all 
1 201 If ) found evidence for strong PAH emission in the 
mid-infrared spectrum of this source and concluded that 
a significant fraction of the bolometric luminosity likely 
results from star-formation as well. Hence, while an AGN 
beamed in the direction of the nebula but obscured from 
our direct view can account for some of the observed dif- 
fuse light, it is not yet clear what fraction of the observed 
Lya and UV continuum emission can be explained by 
this scenario. 

Ultimately, understanding the origin of the diffuse con- 
tinuum components will require better data than are cur- 
rently available. In particular, deep imaging to measure 
the continuum colors, deep polarization observations to 
determine what fraction, if any, of the continuum light is 
scattered, and a map of the velocity field measured from 
a non-resonant line (e.g.. Ha or [OIIIJA5007) are neces- 
sary. Given the apparent faintness and redshift of this 
target, these observations await JWST. 

4.5. Smooth, N on- filamentary Morphology 

Another key morphological result of this work is the 
finding that the diffuse emission components (Lya and 
UV) in LABd05 are remarkably smooth and round. In 
particular, there is no evidence for the kind of bubble- 
like structures that have been ta ken as evidence for the 
superwind outfiow scenario fe.g.. iTaniguchi et al.|[200l[: 
iMori et al.ll2004l : IMatsuda et al.ll2004[ ). The diffuse emis- 
sion is also not particularly clumpy or filamentary, in 
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contrast to the predictions of recent cold flow simula- 
tions that suggest the morphologies of Lya nebulae pow- 
ered by cold accretion should be asy mmetric and n arrow 
with "finger- like extensions" (Gocr dTet al.l[2010l ). In- 
stead, the Lya morphology of LABdOS is quite symmet- 
ric and well-described by an exponential disk. While the 
typical axis ratio of the predicted L ya nebulae appear 
(Figures 7 — 9 of iGoerdt et al.ll2010[) to be in the range 
b/a « 0.25 — 0.5, LABdOS is much less e longa ted with 
an axis ratio of 0.79. The IGoerdt et all (|2010 f) models 
also suggest that clumps associated with the inflowing 
streams should provide an important contribution to the 
total luminosity. Quantitative estimates are not given in 
the paper, but from their Figures 7 — 9 we estimate that 
there should be 3 — 4 significant clumps (with surface 
brightnesses > 0.1 times that of the central peak) within 
the virial radius (w 70 kpc). No significant clumps are 
seen in the Lya emission from LABd05. 

The surface brightness profiles of Lya emission from 
model nebulae are also typically more centrally concen- 
trated than that observed for LABd05. In Figure [9] we 
overplot the predicted surface bri ghtness profile s frorn 
recent cold flow simulations. The IGoerdt et al.l (|2010f ) 
prediction of a r~^'^ power law, which we have scaled to 
match the total observed Lya flux for LABd05 inside a 
radius of 5", is a poor fit to the shape of the observed 
surfac e brightness profile. The iFaucher-Giguere et al.l 
(|2010j) predictions shown are based on their two most 
realistic treatments for Lya emission from gravitational 
cooling (their models #7 and #9, which include pre- 
scriptions for self-shielding). The lower bound of each 
region shown corresponds to the prediction for a fidu- 
cial 2.5 X 10^^ Mq halo mass model at 2: = 3; the upper 
bound of each region is the same profile scaled up to 
a halo mass of 10^'^ M(r) based on the pred icted LLya- 
M/iaio relation (jFaucher-Giguere et al.ll20To() and under 
the naive assumption that the profile shape is constant 
as a function of halo mass. Model #7 can in princi- 
ple reach the peak Lya surface brightness we observe in 
LABd05, assuming a sufficiently massive halo, but the 
profile shape is much more centrally concentrated than 
is observed. Model #9 is orders of magnitude too faint, 
even for the most massive halos. 

We note that although the expected luminosity scal- 
ing is the most basic output from models of Lya neb- 
ulae powered by cold accretion, this has turned out 
to be particularly difficult to predict robustly. Early 
models suggested that the Lya emission from gravita- 
tional cooling sho uld be similar to what is observed i n 
Lya nebulae fe.g.. lYang et al.ll2006l: iGoerdt et al.ll2010( ). 
However, more recent work has argued that these Lya 
nebula luminosity predictions may be orders of magni- 
tude too high due to the effects of self-shielding and 
that cooling radiation alone is an unlikely explanation 
for the Lya emission of the mos t luminous Lya nebu- 
lae (jFaucher-Giguere et al.l [20T0I) . The question of the 
predicted Lya luminosity from cold accretion is still a 
matter of some debate, but interestingly our analysis has 
shown that even if the question of the luminosity scal- 
ing is ignored, key morphological discrepancies remain 
between existing cold flow models and what is seen in 
LABd05. 



4.6. Similar UV and Lya Surface Brightness Profiles 

Due to the effects of resonant scattering, there is a 
generic expectation that high redshift sources of Lya 
emission should be surrounded by low surface bright- 
ness halos of resonantly scattered Lya emission (e.g., 
ILoeb fc Rvbickilll99l fZhens et al.l[2010l) . After a num- 
ber of observational studies uncovered possible hints 
of this extended emiss ion in sample s of Lya-em itting 
galaxies (e.g ., iHavashino et al.l 120041: lOno et all [2010), 
iSteidel et al.l ()2011[ ) used a stacking analysis to demon- 
strate convincingly that extended Lya halos appear to 
exist around all classes of star-forming galaxies. They 
noted the large extent of these Lya halos relative to the 
much more compact UV cores and pointed out the sim- 
ilarities in surface brightness profile shape between the 
Lya halos around stacked continuum-selected LBGs and 
stacked Lya nebulae in the same field. They concluded 
that, if one could image deeply enough {FLya ~ 10~^^ 
erg s~^ cm~^), all continuum-selected LBGs would be 
classified as extended Lya nebulae. The obvious ques- 
tion then becomes: are giant Lya nebulae simply scaled 
up versions of this phenomenon with the large extent 
driven simply by resonant scattering of Lya photons 
from a single or several central sources? 

In the case of LABd05, it does not appear that the 
large Lya extent can be explained simply as a result 
of resonant scattering. The lack of an obvious central 
source of Lya photons is the first challenge (as discussed 
in Section l4.3p . but another inconsistency appears when 
considering the Lya and UV surface brightness profiles 
(Figure [TO)) . Instead of a compact UV core surrounded 
by an extended Lya halo, the UV emission in LABd05 
is nearly as spatially extended as the Lya, a clear indi- 
cation that the the nebula's large size is primarily the 
result of some other mechanism or geometry. Further- 
more, the observed Lya and diffuse UV continuum sur- 
face brightness profiles are remarkably similar in their 
properties. They both show similar radial distributions, 
are well-described by nearly-exponential disks, and have 
approximately elliptical shapes with similar axial ratios 
(Tabled]). In comparison, the stacked UV and Lya pro- 
files derived for continuum-selected star-forming galaxies 
differ significantly (Figure [TOl Stcidcl ct al. 2011). Thus, 
while it seems likely that resonant scattering can explain 
some of the remaining differences between the two com- 
ponents, it does not appear that the large spatial extent 
of LABd05 is solely a result of Lya resonant scattering. 

4.7. A Galaxy Group in Formation? 

LABd05 is unique only because of the existence of 
deep and high-spatial resolution broad- and narrow-band 
imaging data. These data provide many pieces of evi- 
dence in support of the idea that this is a young, forming 
system. There are numerous small, low-luminosity, disky 
galaxies, many of which have very blue colors. Even the 
reddest objects have colors that are consistent with ages 
of less than a few hundred million years, perhaps only 
a hundred million years if they are modestly reddened 
(i.e., E{B — V) fa 0.4). The large and luminous Lya 
halo and the detection of faint, diffuse UV continuum 
emission in the region also suggest that the system is 
energetically young. In addition, the fact that 9 of the 
compact galaxies that are likely to be associated with the 
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system (i.e., roughly half the candidate members, with a 
total luminosity of 1.2 L*) lie within a small projected 
area «30 kpc in diameter suggests that the system may 
also be dynamically young, the dynamical time for this 
region being only ^14 MvrFl 

We speculate that this giant Lya nebula is the progen- 
itor of a galaxy group, witnessed in the process of forma- 
tion. Under this assumption, the luminosity distribution 
of the member galaxies can provide a unique perspective 
on the "initial luminosity function" of galaxies. Spec- 
troscopic redshift measurements of the galaxies within 
the system will be key to confirming membership and 
determining the dynamical mass of the system. While 
we only have spectroscopic redshifts for 3 sources in the 
region, we have argued both from the excess of galaxies 
in the vicinity of the nebula and from the colors of these 
galaxies that most of the compact objects observed are 
likely members of the system. By assuming that all the 
galaxies that lie above the dashed line in Figure |3] are 
members, we have constructed a luminosity function for 
the system, as shown in Figure [14] (bottom panel). Sum- 
ming all the UV luminosity contributed by the candidate 
member galaxies and the diffuse continuum results in a 
total of «23.4 AB mag in the Veoe band, or « 3L* . It 
is therefore possible that this system could evolve into a 
small group, with the smaller galaxies merging into larger 
systems over a few dynamical times. 

Clearly better data are needed, both to confirm the sys- 
tem membership and to measure the stellar masses of the 
member galaxies more robustly. Nevertheless, LABd05 
provides the tantalizing hope that detailed studies of 
more such systems, even statistical studies, can result 
in a determination of the initial mass function of galax- 
ies, analogous to the manner in which studies of stellar 
clusters in our own Galaxy have yielded the stellar ini- 
tial mass function. While LABd05 is only one source, 
we note that there are other Lya nebulae that appear to 
be similar in morphology. For example, one of the Lya 
nebulae in the SSA22 field at z w 3.1 (LABI) shows 
multiple embedded galaxies and a hint of diffuse UV 
continuum in between the galaxies in ground-based data 
(jMatsuda et al. 2 007i ). As discussed in Section l43l the 
Lya nebula found in the GOODS-S field at z ~ 3 . 16 ha s 
no central continuum counterpart (|Nilsson et al.l 120061 ). 
Deep HST imaging of a larger sample of giant Lya neb- 
ulae will be important for understanding the extent to 
which the morphology and galaxy properties observed 
in LABdOS are characteristic of Lya nebula systems in 
general. 

5. CONCLUSIONS 

Using high resolution HST imaging, we have taken a 
census of all the compact sources within a large Lya neb- 
ula at z sa 2.656. We find that the Lya nebula system 
contains numerous compact, young, disky galaxies and 
an obscured AGN that are all located tens of kilopar- 
secs from the peak of the Lya emission and provide a 
negligible contribution to the ionization of the nebula. 
The observed luminosity function shows that the com- 

-"^^ This calculation assumes that L* corresponds to a mass of 
;^2 x 10^° MQ dBh apley et al. 2005; Reddy et al. 2006; Erb ct al] 
[200^; IReddv fc Steidell I200a) and a stellar-to-halo mass ratio at 
this mass of 10-2 for z = 3 HMoster et al.|[20Tcn . 



pact sources within the system are predominantly low lu- 
minosity (^ O.IL*) galaxies, highly suggestive of a galaxy 
forming environment. The large-scale morphology of the 
system is characterized by the lack of a central galaxy 
at or near the peak of the Lya nebula and the presence 
of several nearly coincident, smooth, and spatially ex- 
tended emission components (Lya, Hell, and diffuse UV 
continuum) . These morphological results — in particular 
the lack of a central galaxy and the offset morphology — 
disfavor models of outfiows, cold flows, and resonant scat- 
tering halos, suggesting that while these phenomena may 
be present, they are not sufficient to explain the power- 
ing and the large extent of giant Lya nebulae. Based on 
these observations, we speculate that large Lya nebulae 
are progenitors of low-redshift galaxy groups or low-mass 
clusters. 
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Table 1 

HST Observations of LABdOS 



Instrument 


Filter 


Exposure Time 


Ac 


Bandpass Width 


Restframe A 






(min) 






at 2 Pi 2.656 


HST/ACS 


FR462N ([O ii] outer ramp) 


216 


4448 A'' 


89 A 


LyaA1216 


HST/ACS 


F606W (Broad-band V) 


129 


5907 A 


2342 A 


1295-1936 A 


HST/ACS 


FR601N ([O III] outer ramp) 


129 


5998 A^ 


120 A 


He nA1640 


HST/NICMOS NIC2 


FllOW (Broad-band J) 


281 


1.1 /im 


0.6 /im 


2188-3829 A 


HST/NICMOS NIC2 


F160W (Broad-band H) 


281 


1.6 /im 


0.4 /im 


3829-4923 A 



^ Ramp filters FR462N and FR601N were centered on Lya and HeiiA1640, respectively, at z^2.7 during these observations. 



Table 2 

Astrometric Uncertainty" 



Band Noij^ (Ja cr^ 

(arcsec) (arcsec) 



NDWFS Bw 








Subaru /A445 


5062 


0.10 


0.09 


ACS Lya 


8 


0.17 


0.13 


ACS Hen 


18 


0.12 


0.10 


ACS ^606 


402 


0.12 


0.09 


NICMOS Jiio 


6 


0.12 


0.09 


NICMOS iJieo 


7 


0.12 


0.09 



^ Astrometric uncertainty relative to the NDWFS 

Bw image. 

^ Number of common sources used to compute as- 
trometric correction. 
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Table 3 

Properties of Compact Sources within LABd05 



ID 


Right Ascension 


Declination 


Re'' 




my 


b 




^110 — ^^leo'' 


System'^ 




(hours) 


(degrees) 


(arcsec) 




(AB) 


(AB) 


(AB) 


Member? 


36 


14:34:10.981 


33:17:32.48 


0.39 


19.9 


28.24± 


0.38 


1.62± 0.43 


1.69± 0.22 


Ml (AGN) 


26 


14:34:11.036 


33:17:34.47 


0.17 


1.0 


25.36± 


0.03 


0.09± 0.07 


0.32± 0.10 


Ml (LBG) 


59 


14:34:11.041 


33:17:34.14 


0.27 


1.4 


26.15± 


0.06 


0.03± 0.15 


0.36± 0.20 


Ml (LBG) 


28 


14:34:10.913 


33:17:33.80 


0.15 


2.0 


26.86± 


0.11 


-0.34± 0.37 


l.OOit 0.42 


Ml 


32 


14:34:10.990 


33:17:32.89 


0.31 


1.4 


27.17± 


0.15 


0.20± 0.32 


l.lOit 0.33 


Ml 


33 


14:34:11.025 


33:17:32.70 


- 


- 


27.61± 


0.21 


0.95± 0.30 


l.lOit 0.25 


Ml 


34 


14:34:11.085 


33:17:33.02 


- 


- 


27.98± 


0.30 


1.35± 0.37 


1.30± 0.24 


Ml 


35 


14:34:11.059 


33:17:32.76 


0.29 


2.8 


26.46± 


0.08 


0.72± 0.12 


0.94± 0.12 


Ml 


37 


14:34:11.002 


33:17:32.35 


0.27 


2.2 


26.80± 


0.11 


l.lOib 0.14 


1.20± 0.11 


Ml 


1 


14:34:11.036 


33:17:25.76 


0.31 


2.0 


26.63± 


0.09 


0.96± 0.13 


0.79± 0.11 


M2 


18 


14:34:10.764 


33:17:36.75 


0.15 


0.4 


27.29± 


0.16 


0.67± 0.26 


0.51± 0.30 


M2 


21 


14:34:10.956 


33:17:36.06 


0.26 


0.9 


26.73± 


0.10 


0.68± 0.16 


0.48± 0.18 


M2 


24 


14:34:10.926 


33:17:34.81 


0.23 


0.8 


26.46± 


0.08 


0.34± 0.15 


0.32± 0.21 


M2 


29 


14:34:10.743 


33:17:33.85 


0.09 


1.1 


26.81± 


0.11 


0.31± 0.21 


0.63± 0.25 


M2 


45 


14:34:11.258 


33:17:30.24 


0.10 


0.9 


27.59it 


0.21 


1.71± 0.24 


1.77± 0.11 


M2 


52 


14:34:10.845 


33:17:27.10 


0.11 


0.6 


26.69± 


0.10 


0.24± 0.20 


0.30± 0.29 


M2 


57 


14:34:11.039 


33:17:25.56 


0.08 


0.5 


26.61± 


0.09 


0.94± 0.13 


0.87± 0.11 


M2 


o 
z 


1 /I . Q/1 . 1 n on "7 


oo.l ( .zo.yf5 


nil 


U.4 


28.04± 


0.32 








27 


14:34:10.897 


33:17:34.12 


0.21 


0.2 


27.97± 


0.30 








38 


14:34:10.902 


33:17:32.29 


0.20 


2.6 


27.62± 


0.22 








43 


14:34:11.015 


33:17:31.17 


0.16 


0.2 


27.54± 


0.20 








47 


14:34:10.544 


33:17:29.10 


0.17 


0.7 


27.77± 


0.25 








48 


14:34:10.551 


33:17:28.61 


0.16 


0.8 


27.76± 


0.25 








49 


14:34:11.412 


33:17:29.15 


0.08 


0.7 


27.47± 


0.19 








9 


14:34:10.955 


33:17:38.49 


0.17 


0.9 


27.19± 


0.15 


0.82± 0.22 


O.Olit 0.32 


NM 


13 


14:34:11.251 


33:17:37.90 


0.09 


0.9 


27.18± 


0.15 


0.67± 0.24 


-0.29± 0.45 


NM 


40 


14:34:10.555 


33:17:31.92 


0.11 


0.7 


27.46± 


0.19 


0.79± 0.29 


-0.42± 0.57 


NM 


58 


14:34:10.936 


33:17:36.24 


0.50 


4.1 


27.53± 


0.20 


1.13± 0.26 


0.05± 0.32 


NM 


46 


14:34:10.850 


33:17:29.85 


0.11 


0.9 


25.49± 


0.04 


-O.lOit 0.09 


0.30± 0.13 


NM (2 = 3.2)=^ 


60 


14:34:10.854 


33:17:30.05 


0.13 


1.0 


25.97± 


0.05 


0.05± 0.12 


0.24± 0.18 


NM (z = 3.2)^* 



^ Morphological measurements i?e and n denote the effective radius and Sersic index as measured by GALFIT ("Scction fmi t. 

^ Aperture magnitudes for the Vqqq, Jho, and Hieo bands were computed using 0.4" diameter apertures and aperture corrections 

of [1.15, 1.52, 1.76]. 

^ Membership categories based on optical/NIR colors and proximity (see Section I3.1.3| l: "Ml" - likely system member inside a 
radius of 2.l" from the AGN (#36); "M2" - likely system member but beyond a radius of 2.l" from the AGN (#36); "NM" - likely 
non-member. Objects with membership confirmed by a spectroscopic rcdshift are indicated with parenthetical remarks. Sources with 
no membership designation were not well-detected in one or more bands. 

^ VeoG aperture magnitude uncorrected for contamination by Lya emission at 2 3.2 {Fi^ya — 5.15 X 10~^^ erg s""'^ cm~^; Paper I). 
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Table 4 

Properties of Diffuse Components witfiin LABdOS 



Band 

Restframe wavelengtfi (A) 


f298-f939A 


Diffuse Continuum 
-'no 
2298-3720A 


3829-4923A 


DifEuse Lya 

7A445 
1640-1716A 


Diffuse Lie 11 


GALFIT Parameters'' 


Right Ascension (flours) 
Declination (degrees) 
Size (arcsec) 
Sersic Index (n) 
Axis Ratio (b/a) 
Position Angle (°) 
Total Magnitude (AB) 
Total Luminosity 
(fO^S erg s-l_H"z-i) 
(lO'*^ erg s-i) 


14;34:10.940 
+33:f7:29.87 
Re = 1.67 ±0.08 
0.77 ±0.05 
0.77 ±0.03 
-43.49 ± 5.48 
23.81 ±0.04 

= 17.13 ±0.63 






14:34:10.986 
+33:17:30.43 
Re = 2.48 ±0.03 
0.73 ±0.01 
0.79 ±0.01 
-15.51 ± 1.24 
21.99 ±0.01 

LLy„ = 101.61 ±0.94 


14:34:10.925*= 
+33:17:29.92'' 
FWHM = 0.58 - 1.0^= 



Aperture Photometry'^ 



Magnitude (AB) 


25.39 + 0.09 


> 25.75 


> 25.79 


24.60 ±0.01 




Luminosity 












(10^8 erg s-^Hz-^) 


= 4.00 ±0.35 


< 2.87 


L„ < 2.77 






(10*2 erg s-i) 








LLyc, =9.15 ±0.11 


LHcii = 2.35 ±0.02° 



^ Errors on morphological parameters are formal fitting errors reported by GALFIT but do not include errors resulting from the continuum 
subtraction or previous fitting of embedded compact sources. They therefore likely underestimate the true uncertainty in fitting this complex 
system. 

^ Approximate position of the Hell source measured from ground-based spectroscopic data (Paper I). 

^ Limits on the size derived from ACS Hell imaging and ground-based spectroscopic data, as described in Section |3.2.2| 

^ Aperture photometry measured at the center of the nebula within a small 1.6" diameter aperture, chosen to avoid nearby compact sources. 
No aperture corrections were applied. Limits are 3tT values. 

^ Hell luminosity measured within a 4.5 X 1.5" spectroscopic slit assuming z — 2.6562 (Paper I). 
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Figure 1. Multi- wavelength imaging of LABdOS. In eacii panel, [0",0"] is centered on the location of tlie obscured AGN (see Section l3.1.2t . 
and the position of the Lyo centroid measured from the SuSaru imaging is shown (red cross). The compact components identified in ground- 
based imaging (Galaxies "A" and "B" from Paper I) are labeled on the NDWFS B\y image. The lower right-hand panel shows a zoomed-in 
version of the NICMOS Hiep image with the position of the M IPS 24/im and IRAC 3.6/im source positions indicated (blue and green la- 
error circles, respectively: iDev et al.ll2005l ; [Horiian et al.ll2003 'l. At z ^ 2.656, 1" corresponds to a physical scale of 7.96 kpc. 
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Figure 2. A composite image of LABdOS made by stacking the VboBi ^liOi and Hiao imaging. All compact sources detected above the So- 
limiting magnitude in the Veos band and located within 7" of the obscured AGN (#36, located at [0",0"]) are labeled with the ID number 
used in TablefS] The contours correspond to Lya surface brightness levels of [1, 3, 5, 7, 9]xl0~^'^ erg cm~^ s~^ arcsec"^, as measured 



from the continuum-subtracted Subaru /A445 (Lya) imaging 
shown (red cross). 



The position of the Lya centroid measured from the Subaru imaging is 
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Figure 3. Jiio — ^^160 vs. Vgoe — >^lio color-color plot for compact sources in the vicinity of LABdOS. Known spectroscopic members are 
denoted with large green circles (the LBG system, #26+59, and the counterpart to the obscured AGN, #36). The known spectroscopic 
interlopers (#46+60) are shown as measured (large open circles) and after correcting for the contribution of Lya emission to the Veoe-band 
(open triangles; see Table IS}. The greyscale contours represents all galaxies from the HUDF above the magnitude limits of our V"606i JiiOi 
and J/i60 data; the line contours represent the subset with photometric redshifts at the redshift of LABdOS (zphot = 2.656 + 0. 15; ICoe et al.l 
1200611 . The dashed black line indicates the division used for membership assignment (Section 13. 1.31 1. The appropriate reddening vector is 
shown for E(B — V) = 0.1 computed at 2: fs 2.656. 
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Figure 4. Jiio — Hiqq vs. Veoe — Jiio color-color plot for compact sources in the vicinity of LABd05. Member sources in category 
"Ml" are shown as large filled circles — with known spectroscopic members in green, i.e., the LBG system (#26+59) and the obscured 
AGN (#36) — and those in "M2" as small filled circles. Non-member sources ( "NM" ) are shown as open circles; the known spectroscopic 
interlopers (#46+60) are shown as measured (large open circles) and after correcting for the contribution of Lya emission to the Veoe-band 
(open triangles; see Table (3)1 . A scries of single stellar population model age tracks (unreddened burst, solar metallicity, spanning burst 
ages of 5-1400 Myr; FEruzual iSTCha rlot 2003) are overplotted for different redshifts. The thick green line corresponds to z ^ 2.656, the 
redshift of LABd05. The 25 Myr, 100 Myr, and 320 Myr positions along the tracks are indicated with small filled circles. For comparison, a 
constant star-forming model (dotted green line) and a low metallicity model (Z = 0.0001, dot-dashed green line) are shown, both sp annin g 
the same age range at the redshift of LABd05. The dashed black line indicates the division used for membership assignment (Section l3. 1.31 1. 
The appropriate reddening vector is shown for E{B — V) = 0.1 computed at 2 2.656. 
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Figure 5. Photometric redstiift distribution for galaxies in the HUDF that satisfy the color cut used for membership a ssignment in thi s 
work (black solid histogram) in comparison with the distribution for the full HUDF galaxy sample (blue dotted histogram: ICoe et al.l l2006'). 
This comparison demonstrates that the proposed color cut is successful at selecting high redshift sources (90% at z > 1) with the peak of 
the resulting redshift distribution centered on the redshift of LABdOS (red dashed line). 
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Figure 6. LABd05 galaxy membership. The obscured AGN is centered at [0",0"] in both panels, and the position of the Ly« centroid 
measured from the Subaru imaging is shown as a red cross. Left: The ACS Veoe image, with circles representing sources that have been 
flagged as members of the system based on optical/NIR colors and dis tance from the AGN ("Ml" spectroscopic members: filled green 
circles, "Ml": large black circles, "M2": small black circles; see Section l3.1.3| l. Right: The continuum-subtracted Subaru IAAA5 (Lyo) 



imaging, with contours at Lya surface brightness levels of [1, 3, 5, 
offset by >1.9"~15 projected kpc from all of the member galaxies. 



7, 9]xl0 erg cm s arcsec . Note that the Lya emission is 
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Figure 7. GALFIT parametric fit to the Vgoe imaging of LABdOS. Top: The original Veoe image and GALFIT model. Bottom: 
The residual image and smoothed residual image (FWHM = 10 pix = 0.5" Gaussian kernel) reveal a diffuse UV continuum component. 
Note that the diffuse Vgoe continum emission is nearly coincident with the Lya centroid (red cross; an offset of ^ 0.8"), and that both 
components are in turn offset by fsl.9-2.6"~15-21 projected kpc from the position of the obscured AGN at [0",0"]. 
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Figure 8. Lya imaging of LABd05. The obscured AGN is centered at [0",0"] in all panels, and the position of the LyQ centroid measured 
from the Subaru imaging is shown as a red cross. Left: The continuum-subtracted Subaru /A445 (Lya) imaging, with contours at Lyo 
surface brightness levels of [1, 3, 5, 7, 9]xl0~^^ erg cm~^ s~^ arcsec"^. Middle: The ACS continuum-subtracted Lya image. Right: 
The ACS continuum-subtracted Lyo image smoothed to match the PSF of the ground-based Subaru imaging (FWHM=0.7"). No compact 
knots or clumps are detected in the ACS Lya imaging. 
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Figure 9. Surface brightness profiles of the Lya and diffuse UV continuum emission in LABd05. Top: Lya surface brightness profiles 
computed using the Subaru Lyo (filled circles) and ACS Lya (small open circles) imaging as well as the GALFIT parametric fit to the 
Subaru Lya imagin g (red line). P redictions from cold flow simulations are overplotted and labeled with 'GIO' (dashed purple line; the 
prediction from Goe rdt et al.ll20ldl . sc aled to match the total Lya luminosity of LABdOS within a 5" radius) and 'FGIO' (hatched blue 
regions bounded by dot-dashed lines; [Faucher-Giguere et al.|[2010h . as discussed in Section 14.51 The observed profile of LABd05 is not 
well-described by the predictions from existing cold flow simulations. Bottom: Veoe surface brightness profile of the diffuse UV continuum 
emission in LABdOS (squares) as well as the GALFIT parametric fit (green line). A GALFIT parametric fit to the Subaru Lya imaging 
from the top panel is reproduced for reference (dotted red line), after being scaled to match the total observed flux of the diffuse UV 
continuum within a radius of 3". The centroids of the Lya and diffuse UV continuum emission are offset by 0.8". 
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Figure 10. Surface brightness profiles of the Lyo and diffuse UV continuum emission in LABd05, as shown in Figure [9] Top: The 
stacked Lyo surface brightness profiles of continuum-selected galaxies ('LBG-Sll'; blue line) and of Lyo nebulae (' LAB-Sll'; orange li ne) 
are overplotted. Bottom: The stacked continuum surface brightness profile of continuum-selected galaxies from ISteidel et al.l 120111) is 
overplotted {'LBG-Sll'; purple line). The stacked LBGs show compact UV emission surrounded by an extended Lyo halo, likely evidence 
for resonant scattering of Lyo photons generated in the core; in contrast the Lyo and UV surface brightness profiles for LABd05, while 
not identical, are strikingly similar in shape and extent. 
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Figure 11. Limits on the size of the Hell-emitting region derived from the continuum-subtracted ACS Hell imaging. In all panels, the 
position of the Lya centroid measured from the Subaru imaging is shown (red cross), and the approximate position of the Hell emission, 
as measured from ground-based spectroscopy, is indicated with black open crosshairs at [-0.9", -2. 5"]. Top: ACS Veoe and continuum- 
subtracted Hell imaging. Bottom Left: Ratio of the recovered flux of simulated Hell sources with varying sizes (i^simuiated) divided by 
the Icr limiting flux of the image (o-^f-y), shown as a function of the FWHM of the simulated source (0.4" diameter apertures). The red 
band denotes the full range of results from 100 Monte Carlo trials. Fsimulated/'^sky = 3 is shown (dotted line) along with the derived lower 
limit on the si ze of the Hell-emitting region (FWHM > 0.58"; dashed line). Any source smaller than this would have been detected at > 3cr 
(Section l3.2.2l l. Bottom Right: A simulated Hen image containing a representative set of model sources (FWHM = [0.1", 0.2", 0.4", 0.6"]). 
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Figure 12. Number counts for sources detected above the Scr limiting magnitude in the Vgoe band that are located within a 7" radius of 
the AGN in LABdOS (black solid line) and over the entire 207" X 205" ACS pointing (dashed black line). The 80% and 50% completeness 
limits are indicated (thin and thick dotted lines, respectively). The Vfjos-band number counts from the HD F-N a nd HDF-S (red and green 
lines; WiUiams ct al. 1996; Casertano et al. 2000), the HUDF (purple line; Bcckwith et al. 2006; Coc et aL|[200i), and the GOODS-N and 
GOODS-S fields (light and dark blue lines: iGiavalisco et al.i,2004 ) are shown, remeasured using aperture magnitudes consistent with our 
analysis. 




Figure 13. Overdensity scale factor, relative to the field number counts from HUDF, of galaxies within a 7" radius of the AGN that 
are detected above the 5cr limiting magnitude in the Veoe band (black solid line). The LABd05 region is overdense by at least a factor of 
^ 4 relative to the field. The 80% and 50% completeness limits (thin and thick dotted lines, respectively) and the result of ap plying the 
completeness correction are shown (red open circles). The top axis is labeled in terms of L* at z « 3 (A/* = — 20.8: rReddv et~al.ll200SV 
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Figure 14. Luminosity function of galaxies within LABdOS. Top: Sources detected above the 5cr hmiting magnitude in the Veoe band 
and located within a 7" radius of the AGN (black open histogram) and corrected for incompleteness (thin dashed histogram). The 80% and 
50% completeness limits are shown (thin and thick vertical dotted lines, respectively). Bottom: Sources detected above the 5<j limiting 
magnitude in all three bands and located within a 7" radius of the A GN ( black open histogram). The filled and hatched histograms 
represent galaxies in groupings "Ml" and "M2", respe ctively (Section 13. 1.31 1. The thick dashed curve in both panels is a "statistical" 
lumin osity function for LABdOS, as described in Section l4.2.2l The top axis is labeled in terms of L* at 2: fs 3 (Af* = — 20.8: iReddv et al.l 
[2003) . 
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Figure 15. Veoe morphologies derived using GALFIT for galaxies near LAB dOS. O pen histograms contain all sources; solid and hatched 
histograms represent sources in groups "Ml" and "M2", respectively (Section 13. 1.31 1. Top: The histogram of galaxy effective radii (Re), 
with the effective radius of the ACS PSF shown for reference (dotted line). Bottom: The histogram of Sersic indices (n), with the values 
for exponential (n = 1) and De Vaucouleurs (n = 4) profiles shown for reference (dotted lines). 
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APPENDIX 
POSTAGE STAMPS 

Postage stamp images of the compact sources in the vicinity of LABdOS are shown in Figures [TBll?!] From left to 
right, the panels display the stacked image and the individual Veoe, Jiio^ and iJigo-band images. Object ID numbers 
are shown in the upper left-hand corner of the stacked image, and the membership category (where applicable) is given 
in the upper right-hand corner. The AGN (#36) and both components of the LBG (#26 and #59) are labeled in the 
lower right-hand corner. 

We note that the current data are not sufficient to distinguish whether the object pairs (#26+59, #1+57, #21+58, 
and #46+60) are true companions or just morphological peculiarities (i.e., tidal features, dust lanes, etc.) associated 
with the primary object. For our analysis we have chosen to treat each of the object pairs as two separate objects. 
Treating them each as a single object would decrease the number of sources in the "Ml" category by one and the "M2" 
category by two, but this does not significantly alter our conclusions regarding the nature of the member galaxies or 
the overall luminosity function within the LABd05 system. 
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Figure 16. Postage stamp images of the compact sources in the vicinity of LABd05. 
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Figure 18. Postage stamp images of the compact sources in the vicinity of LABd05. 




Figure 19. Postage stamp images of the compact sources in the vicinity of LABd05. 
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Figure 20. Postage stamp images of the compact sources in the vicinity of LABd05. 
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Figure 21. Postage stamp images of the compact sources in the vicinity of LABd05. 



